OUTP-04/17P 



Asymmetric Sneutrino Dark Matter 
and the £\/^dm Puzzle 



Dan Hooper 1 , John March- Russell 2 , and Stephen M. West' 

1 Astrophysics Dept., University of Oxford, 
Denys Wilkinson Building, Oxford 0X1 3RH, UK 

2 Rudolf Peierls Centre for Theoretical Physics, 
University of Oxford, 1 Keble Road, Oxford 0X1 3NP, UK 

(October, 2004) 



Abstract 

The inferred values of the cosmological baryon and dark matter densities are strikingly 
similar, but in most theories of the early universe there is no true explanation of this fact; 
in particular, the baryon asymmetry and thus density depends upon unknown, and a pri- 
ori small, CP-violating phases which are independent of all parameters determining the 
dark matter density. We consider models of dark matter possessing a particle-antiparticle 
asymmetry where this asymmetry determines both the baryon asymmetry and strongly 
effects the dark matter density, thus naturally linking £\ and fW- We show that sneutri- 
nos can play the role of such dark matter in a previously studied variant of the MSSM in 
which the light neutrino masses result from higher-dimensional supersymmetry-breaking 
terms. 



1 Introduction 



For some time it has been apparent that the inferred values of the cosmological baryon 
and dark matter densities are strikingly similar. For example, the latest WMAP- 
determined range for the dark matter density, 0.129 > ^cdm^ 2 > 0.095, is within a 
factor of a few of the combined WMAP and big-bang nucleosynthesis determined value 
of the baryon density, 0.025 > Q^h 2 > 0.012 P3I2]- in the vast majority of models of the 
early universe, the cosmological baryon and dark matter densities are independently de- 
termined. The surviving baryon density is set by a baryon asymmetry generated during 
baryogenesis, and thus depends upon unknown, and a priori small CP- violating phases, 
as well as unknown baryon- number violating dynamics. In contrast, the dark matter 
density is set by the 'freeze-out' of the interactions that keep the dark matter in equi- 
librium, and is independent of the dynamics of baryogenesis. In particular, although 
weakly-interacting massive particles at the TeV-scale naturally have a relic density of 
(9(1) times the critical density, this is not at all the case for baryons. The Boltzmann 
equations show the size of the baryon relic density in the absence of an asymmetry is 
(9(10 -11 ). Thus the comparative closeness of the baryon and dark matter densities poses 
a puzzle. One possible approach to this problem is to more closely integrate the dynam- 
ics of baryogenesis with that of the origin of dark matter. In particular, it is natural to 
consider models where the dark matter and baryon sectors share a quantum number, 
either continuous or discrete, which provides a relation between their surviving number 
densities and thus energy densities. 1 

Specifically, in this letter we consider models of dark matter possessing a particle- 
antiparticle asymmetry where this asymmetry strongly effects the dark matter density, 
and through the electroweak (EW) anomaly, determines the baryon asymmetry, thus 
naturally linking and ^cdm- (A noteworthy early attempt along these lines which 
shares some features with our model is contained in Ref. [3].) The shared quantum 
number clearly requires that the dark matter particle not be it's own antiparticle, so the 
putative dark matter candidate cannot be the standard LSP neutralino of the broken 
Minimal Supersymmetric Standard Model (MSSM). In fact, we will show in sections 3 
and 4 that sneutrinos can play the role of such dark matter in a previously studied 
variant of the MSSM in which the light neutrino masses result from higher- dimensional 
supersymmetry-breaking terms [UJ |Hl El UHj • This model preserves all the successes of 
the MSSM, such as stability of the weak scale and unification of gauge couplings, while 
being, at least in part, testable at the LHC. 

Before we proceed to the calculation of the relic dark matter and baryon densities 
and the details of our model, it is useful to consider some aspects of the idea of a shared 
quantum number determining the ratio of baryon to dark matter. First, consider the 
simplified case in which, to a very good approximation, the dark matter and baryon 
sectors cannot exchange their conserved quantum number after its first production. As 
an example, consider the situation in which the asymptotic baryon number B = 1 states 
of the Standard Model (SM) have global charge, q, while the lightest asymptotic states 

*An alternate approach is to invoke a form of the Anthropic Principle, see, for example, Ref. 
Another recent idea is explored in Ref. 0]. 
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(with mass m^m) in a hidden sector carry charge, Q. Then conservation of global charge 
implies that q(n h — n^) = —Q(rid m — n-^), where the n's are the number densities 
of the indicated (anti)particles. Further assume that interactions in both sectors are 
strong enough such that almost all antiparticles are eliminated by annihilation with 
their particles (assuming Q/q < for simplicity). This implies 

n b = cn dm with c=\Q/q\, (1) 

which in turn leads to 

7j— = = c • (2) 

Thus the energy densities are indeed related, but the ratio is only naturally 0(1) if the 
ratio of the baryon to dark matter masses is not too small. It is important to note that 
this disfavours models of the above type where the dark matter candidate is a 'hidden 
sector' particle and favours a particular class of models where the dark matter candidate 
has weak scale mass, its mass arising either from electroweak symmetry breaking, or from 
the dynamics which drives electroweak breaking, such as supersymmetry softly broken 
at the weak-scale. 2 

As a consequence, the naive statement of Eq.(^) also requires modification as the 
approximation of negligible interactions between the two sectors, and thus negligible ex- 
change of the shared conserved quantum number, does not hold in realistic models. In 
this situation, the ratio of the conserved quantum number stored in the two sectors, c, 
is not simply given by Q/q, but is determined by the 'chemical' equilibrium conditions 
between the two sectors just before the freeze-out of the relevant interactions. We will 
further discuss this point in section 2. Moreover, two other alterations to the naive re- 
lation Eq.(J2J) can be present. First, the shared quantum number may only be conserved 
modulo n, in which case 'self annihilations may occur, and second, it is not always the 
case that dark matter interactions are efficient enough such that n-^ <C fidm- Neverthe- 
less, as we argue in the next section, a careful consideration of such models shows that 
a version of the naive statement Eq.(J2J) can indeed hold. 



2 The Effect of a Matter- Antimatter Asymmetry on 
the Density of a Thermal Relic 

If there were an asymmetry between the density of dark matter particles and dark anti- 
matter particles in the early universe, it may have a considerable effect on the relic 
density calculation [TT]. In this section, we review the relic density calculation for a 
single species, two species with no asymmetry and two species with an asymmetry. 

If dark matter consists of a single species, the relic density calculation is straight- 
forward. The most common example of a single species dark matter candidate would be 

2 In extensions of the SM such as the MSSM where many gauge-non-neutral fermionic and scalar 
states also gain mass at this scale, the beta-function coefficient for SU(3)qct> sharply increases as 
energies are reduced below the weak scale, and the confinement scale and thus baryon mass is naturally 
not very much smaller than m wea k. 
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a neutralino LSP with negligible coannihilations with other SUSY particles. In such a 
case, the resulting relic density is simply 



, 1.04 x 10 9 x F 



VLh 2 ^ — ~ r, (3) 

Mpi ^gZ(a + 3b/x F y v ; 

where M-p\ is the Planck mass, xp is the WIMP's mass over the freeze-out temperature, 
g* is number of relativistic degrees of freedom available at freeze-out and a and b are 
related to the WIMP's self-annihilation cross section by ov = a + bv 2 + 0(v 4 ), v being 
the relative velocity of the annihilating WIMPs. For a weakly interacting relic, xp is 
generally about 20 and varies only mildly with changes in the cross section or mass. 

If we instead consider a dark matter candidate and its antiparticle, the calculation 
becomes somewhat different. Labeling these species 1 and 2, we now have each particle's 
self-annihilation cross section, an = a~22, and the particle-antiparticle annihilation cross 
section, 012 = 021 ■ In the limit that o\i is much smaller than an, we simply get two 
decoupled species which each provide the relic density of Eq.©. Otherwise, the effective 
cross section becomes an + 012 for each species. For the two species together, a density 
of 

„ t2 n 1.04 x 10 9 x F 1 

Qh 2 ^2 x - (4) 

Mpi v ^r(a 1 i + ai a + 3(6u + 6i 2 )/zf) 

is produced. This case is essentially the calculation of the density of a relic with the 
coannihilation of a degenerate particle [T2] . 

If an asymmetry exists between the number density of particles and antiparticles in 
the early universe, again this calculation becomes modified. The most well known exam- 
ple of this is the baryon-antibaryon asymmetry which leads to the current baryon density 
of the universe despite the very large baryon-antibaryon annihilation cross section. In 
such a case, where the particle-particle and antiparitcle-antiparticle annihilation cross 
sections (011 and 022) are zero, or negligibly small, the minimum relic density can be 
related to the magnitude of the particle-antiparticle asymmetry: 

= ^Lv^ — — , (5) 



^bary 



where Qh1 &j:y is the current density of baryons in the universe, A and Ab ary are the 
particle-antiparticle asymmetries of our relic and of baryons, defined by A = (n — n)/n. 
^bary is known to be of the order of 10~ 9 . m and mb ary are the masses of our relic and 
of baryons (i.e. the proton mass). Of course, if the particle-particle annihilation cross 
section for the relic is not negligible, Eq. (0) will not hold. There will be a generic tendency 
for the density of a relic to move towards this value as a result of an asymmetry, however. 

To accurately assess the effect of an asymmetry on the thermal relic abundance, the 
effective cross sections must be determined by integrating over the thermal history. In 
particular, for each species separately, the densities predicted are 

Q ^ 2 ^ 1.04 x 10 9 x F 1 

Mpi Jgl (aii >e ff + ai2, e ff + 3(&u ieff + bi2 jeS )/x F ) 
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and 




„ , o 1.04 x 10 9 x F 1 

^ 2 = 77 F=7 ; 7^777 -r r~/ r > (7) 



-^Pl \/#* (a22,eff + a21,eff + 3(6 2 2,eff + b 2 l,^[)/xp 

where the effective cross sections are given by 

a ll,eff = All j a 12,eff = %fJ^^T 

022,eff = a 22, a 21,eff = XpJ^^T 

&ll,eff = &11> &12,eff — l 

W = & 22 , &2Mf = 24/~^^±4Jete. 

The factor, e~ x /(e~ x + A), is the depletion of the second species due to the asymmetry 
at a temperature, T = m/x. The effect of this depletion depends on the size of the 
asymmetry and the cross sections involved. If the estimate of Eq. © is considerably less 
than the standard two-species result of Eq.(J3J), then the asymmetry has little impact 
and the sum of Eqs.(JB} and (jZj) is nearly equal to Eq.(j3J. If, however, the estimate of 
Eq.fjSJ) is larger than the standard two-species result, the relic density moves from the 
value found with Eq.Q towards that of Eq.@. If crii/cr 12 is small, then the result will 
be closer to Eq.fjSJ). If (J\\l<J\i is large, then the density will be less modified. Note that 
in some cases, the effect of an asymmetry can be to lower the density of a relic by up to 
a factor of two. This occurs when o\\ and an are comparable thus efficiently deplete the 
second species while leaving the first largely uneffected. 

Although these effects are rather model dependent, the net effect can be large regions 
of parameter space which predict the observed density of dark matter for a WIMP with 
a mass near m ~ m p fiwMAph 2 /flbar y h 2 , for a matter-antimatter asymmetry similar to 
that for baryons. For concreteness, we consider a specific example of such a model in the 
following section. 



3 An Example: Mixed Sneutrino Dark Matter 

Within the context of the Minimal Supersymmetric Standard Model (MSSM), sneutrinos 
do not make a very appealing dark matter candidate. They annihilate quite efficiently, 
resulting in a relic density smaller than the observed dark matter density for sneutrinos 
lighter than about 500 GeV. Furthermore, their elastic scattering cross section is suffi- 
ciently large to be easily observed by direct dark matter experiments. Thus it is natural 
that in contrast to other supersymmetric dark matter candidates, sneutrino dark matter 
has received considerably less attention [TH] . 

Looking beyond the vanilla MSSM, however, it is quite easy to construct a viable 
model for sneutrino dark matter. For example, sneutrinos could mix through /L-terms 
with a singlet forming a mixed sneutrino state. This would reduce the annihilation cross 
section, potentially providing the appropriate quantity of dark matter, as well as reduce 
the elastic scattering cross section to experimentally acceptable levels. In fact, such a 
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scenario has previously been discussed as a way to reconcile the results of the direct dark 
matter detection experiments DAMA and CDMS [8 j . These models have light sneutrino 
mass eigenstates which are mixtures of left-handed 'active' sneutrinos, and right-handed 
'sterile' sneutrinos. 

Such mixing naturally arises in models in which the light neutrino masses are gen- 
erated by supersymmetry breaking effects [3 13 El ED] • The importance to us of these 
models is that the light sneutrino states share a non-anomalous (B — L)-symmetry with 
the baryons which is only weakly broken by the Majorana neutrino masses. In analogy to 
the Giudice-Masiero mechanism for the Higgs //-term there exists a global symmetry that 
forbids a direct mass, MrNN, and neutrino Yukawa coupling, XLNH U , involving the SM 
singlet right-handed neutrino superfield N, but allows non-renormalizable couplings to 
a supersymmetry-breaking hidden-sector field T. The right-handed neutrino mass and 
the neutrino Yukawa coupling, XLNH U , then only arise from l/M planck -suppressed terms 
involving the fundamental supersymmetry breaking scale mi ~ [vM p i) 1 / 2 , where v is the 
weak scale. The relevant terms in the Lagrangian are 

/• / Tt _ T^T T 7 tT"T't \ r / T \ 

c = / dH [ h M NN + h jp N ' N + hB n^ NN +■■■)+ J d2 ° \^m lnh -) ' ( 9 ) 

where all dimensionless couplings g, h, etc, are taken to be 0(1). The structure of the 
Lagrangian can be enforced by some simple symmetries [HI 011]], including R-parity, which 
leads to a stable lightest supersymmetric particle (LSP) dark matter candidate. After 
supersymmetry is broken in the hidden sector at the scale mi, the F and A components of 
the field T acquire intermediate scale vevs, (Ft) ~ (At) 2 ~ irij- Substituting these vevs 
into the Lagrangian shows that after SUSY breaking, there arises: (1) a neutrino Yukawa 
coupling of size ~ 1CT 7 — 1CT 8 , (2) right-handed neutrino and sneutrino masses at the weak 
scale, (3) a weak-scale trilinear scalar A-term, and (4) a right-handed sneutrino lepton- 
number violating 5-term with mag nitude B 2 ~ (few x 100 MeV) 2 . These together lead 
to two sources of neutrino masses, a tree level see-saw contribution, as well as a dominant 
1-loop contribution, leading to a consistent and attractive neutrino phenomenology jHHSJ. 

More relevant for us, the sneutrino mass matrix squared is 

( * - ) ( Kt" )(:•)• (10) 

where we have ignored the small lepton number violating I?-term for the moment. There- 
fore mixing between the left-handed 'active' sneutrino, v, and the right-handed 'ster- 
ile' sneutrino, ft, occurs, producing the light mass eigenstate V\ = —usinO + n*cos8, 
where tan 29 = —\f2Avsm. (5 / (m 2 R — m 2 L ), with m\ = Mjj + m 2 Rsoit and m\ = 
m\ cos 2/3/2 + m 2 Lsoft . Since the coupling of the lighter sneutrino eigenstate to the Z 
is suppressed by sin#, the direct LEP experimental constraints are weakened, and the 
calculation of the sneutrino relic abundances is modified. A numerical analysis of the 
susy parameter space shows that a mixing sin 9 ~ 0.3 is quite natural. 3 (The B 2 mass- 
term splits the CP-even and CP-odd components of i>\ by the amount 5 ~ cos 2 9B 2 jm\. 

3 Note that the size of the A terms is limited by considerations of vacuum stability |10| . 
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This can have important consequences for detection, potentially making consistent the 
DAMA and CDMS experiments depending upon parameters |S] . For the purposes of our 
discussion this splitting is irrelevant.) 

We now turn to the calculation of the relative asymmetry in the sneutrino and 
baryon sectors. The method is a simple adaption of the standard 'chemical' equilibriation 
techniques applied to, for example, the calculation of the ratio B/(B — L) in the SM [T3] 
or MSSM jTH] in the presence of anomaly-induced baryon number violating processes in 
the early universe. For each species of particle we introduce a chemical potential, [Aa, 
and for each interaction in equilibrium, say A + B <-> C + D there is a corresponding 
relation among chemical potentials \la + A*b — fJ>c + V-D- Note that when in thermal 
(not necessarily chemical) equilibrium the CP-conjugate species A has 11-^ = —\ia- The 
asymmetry in particle-antiparticle number density for the relevant case of \x <C T and a 
species of mass m is given by n — n = gT 2 fi/6 for T ^> m, and is Boltzmann suppressed 
by exp(— m/T) for T < m. 

The electroweak B + L anomaly interactions provide an efficient mechanism to trans- 
mute between B- and L-number densities, at least above the EW phase transition tem- 
perature T c where they are rapid. Our task is to calculate the ratio of sneutrino to baryon 
asymmetries at T c where the exchange between the two sectors switches off. (We implic- 
itly assume that the EWPT is 1st order so that anomaly-induced interactions drop out 
of equilibrium below T c . For realistic values of T c the results of our calculation are not 
particularly sensitive to this assumption.) Above the EWPT temperature the chemical 
potentials for all SM gauge bosons vanish, so, for example, members of an EW doublet 
have the same \i. Moreover, efficient flavour mixing among quarks induced by EW inter- 
actions such as uu + 0i <-> URji involving the neutral up- like higgs, or du + 4>2 ^Rj-> 
involving the neutral down-like higgs, ensures that only one chemical potential, say that 
for ul, is required to describe the quark sector. Similarly, efficient flavour mixing among 
the light neutrinos holds for mixing angles, 8ij > 10~ 7 , so that only one chemical poten- 
tial, (j,„ L , is needed to describe the neutrino and sectors. In addition non-perturbative 
EW processes can be thought of creating, say, a neutron made of lhd quarks, and a lhd 
neutrino out of the vacuum for each generation. When this process is in equilibrium it 
imposes an additional relation which can be written 

3fi UL + fi UL = 0. (11) 

When the MSSM susy spectrum is included the analysis becomes more involved. To 
simplify the analysis let us first assume at a temperature T (with T > T c ) that the 
MSSM susy spectrum, including k rhd sneutrinos can be considered light (m < T). As 
shown in Ref. ^B] the usual soft-susy-breaking terms and supersymmetric partners of the 
SM interactions lead to enough relations among the chemical potentials that only 3 are 
needed for an independent set: liq = (/^ 10 — /^ 20 )/2 in the neutral Higgs sector, lm Ul , and 
fi UL . In our case we must add an additional k light (m < T c ) rhd neutrinos and sneutrinos 
with their associated potentials \i VRi and fiv m , but we also have A-term interactions 
between the lhd and rhd sneutrinos, and the very small lepton-number violating rhd 
sneutrino B-term. These together enforce the conditions 
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where the first equality arises from A-term mediated interactions (as the A terms are 
0{m w ) it is simple to show that these are in equilibrium), and the second equality arises 
from SU(2) W x U(l)y gaugino interactions. 4 In terms of //o,£Vl, an d Hv L , the condition 
Eq.()ll|l is unchanged in the case of the MSSM + k rhd neutrino supermultiplets. 

The final step in the calculation of the relative number densities proceeds by ex- 
pressing the relevant densities in terms of the independent chemical potentials. For the 
baryon number density B the result is simply 

B = 6T 2 /v- (13) 

For the sneutrino asymmetry the result depends upon the mass spectrum of the sneutrino 
mass eigenstates relative to T c . In the range of parameters of interest to us where the 
light sneutrino eigenstate u% is the LSP we have mi < T c . On the other hand the heavier 
eigenstate 1)2 can have m 2 > T c or m 2 < T c . If both eigenstates satisfy m < T c then the 
sneutrino density is (allowing for a possible k such states) 

kT 2 2kT 2 

S = —(^i>L + »Ur) = —^-VVL ( 14 ) 

while if only mi < T c then S is half as large. Finally using the non-perturbative relation 
Eq.JIU we find tnat 

t (1^) 

^bary B 3 6 

depending on m 2 . We view k = 1 as most likely, so we specialize to that case in the next 
section. Finally we mention that a full susy-mass-threshold dependent calculation of the 
relative asymmetries introduces only small corrections to this result (see Ref.|16j for a 
similar analysis in the MSSM case). 

In the end the vital point is that independent of whatever dynamics at scales E > T c 
produces an asymmetry in either the baryon sector, or in the neutrino or sneutrino 
sector, the (B + L)-anomaly-induced interactions together with EW gaugino and A-term 
interactions automatically distributes the asymmetry between the baryons and the light 
sneutrino states, with a predictable A/A^j ratio. The asymmetry could originate from, 
e.g. a GUT-based baryogenesis mechanism, or maybe more interestingly in the context of 
the sneutrino dark matter model there is the possibility that it could arise via a resonant 
leptogenesis mechanism at the TeV-scale as discussed in Ref.[10j. The end result is that 
we expect 1/3 > A/A^y > 1/6 independent of the source of the asymmetry. 



4 Asymmetric sneutrino and baryon abundances 

We will now apply the formalism of section El to the sneutrino model. Given the value of 
1/3 > A/Abary ^ 1/6, we expect the role of the matter-antimatter asymmetry to be most 

4 Note that the rhd neutrinos only interact via the small neutrino Yukawa, A ~ 10~ 7 — 10~ 8 , so 
these states are not in equilibrium at T ~ T c . However, fi l/R does not enter our final expression for the 
sneutrino to baryon number density ratio, so its value is irrelevant. 
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important for sneutrinos with masses on the order of m ~ few x m p fiwMAph 2 /flbaryh 2 — 
30 GeV. For this reason, we will focus on the mass range of about 10 to 80 GeV. In this 
range, the most important annihilation channel is s-channel Z-exchange to fermions. 
This is a p-wave amplitude, and has no a-term in the expansion in v 2 (see section |2J). 
This is also only a sneutrino-antisneutrino process (i.e. crn^, not cxi2,2i)- Other channels 
include t-channel neutralino and chargino exchanges to neutrinos and charged leptons. 
The channel to charged leptons includes both a and b terms in the v 2 expansion. Addi- 
tionally, sneutrino-sneutrino pairs (or antisneutrino-antisneutrino pairs) can annihilate 
via t-channel neutralino exchange to neutrino pairs of the same flavor. This is the only 
contribution to 012,21 in this model. 

The annihilation cross sections for the neutralino and chargino exchange processes 
depend on the masses and mixings of these exchanged particles, although the Z-exchange 
channels do not. We have adopted the parameters M!=300 GeV, M 2 =300 GeV, /i=600 
GeV and tan (3 = 50 throughout our calculation. Additionally, there is the mixing angle 
between the sneutrino and the singlet which can be varied. This quantity has the effect 
of reducing all electroweak couplings of the sneutrino eigenstate by a factor of sin 2 9, i.e. 
reducing the annihilation cross sections by a factor of sin 4 9. 

Finally, one can consider the s-channel exchange of Higgs bosons to bb (or ZZ and 
W + W~ above the appropriate mass thresholds). This process is not suppressed by sin 4 9. 
The efficiency of this process depends on the magnitude of the trilinear scalar coupling 
and is only important near the Higgs resonance (m„ ~ m^/2) or above the gauge boson 
mass thresholds. 

Our results are shown in figure H The shaded regions of the parameter space predict 
a relic density within the range measured by WMAP (0.129 > ^cdm^ 2 > 0.095). In 
the left frame, no matter-antimatter asymmetry was included. This frame can be easily 
understood as requiring greater suppression of electroweak couplings (smaller sin 9) as the 
sneutrino mass approaches the Z mass. The dip at 56-59 GeV is due to s-channel higgs 
exchange to bb. In the right frame, a matter-antimatter asymmetry of A/A^j ~ 1/6 was 
included. To further illustrate this effect, we show the result of this calculation across 
one value of sin 9 in figure El Below about 30 GeV, the matter-antimatter asymmetry 
has little effect on the calculation and the solid and dot-dashed lines fall nearly on top 
of each other. In the range of roughly 30-70 GeV, however, the asymmetry pulls the relic 
density above the standard symmetric result into the range favored by WMAP. Above 
this range, sneutrino-antisneutrino annihilation decreases, leading to larger relic densities 
for the case with no asymmetry. The relic density for the asymmetric case, however, is 
largely determined by the sneutrino-sneutrino annihilation cross section in this region, 
so does not increase as rapidly, therefore resulting in a relic density much closer to the 
preffered value, even for > 70 GeV. 

5 Indirect Detection 

Sneutrinos in the mass range we have studied here annihilate dominantly through s- 
channel Z exchange in the early universe. This annihilation channel consists of a 6-term 
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Figure 1: The regions of parameter space which provide the quantity of mixed sneutrino 
cold dark matter measured by WMAP, 0.129 > fi C DM^ 2 > 0.095. In the left frame, the 
standard calculation with no matter-antimatter asymmetry is used. In the center frame 
the effect of a matter-antimatter asymmetry with A/A^ ary ~ 1/6 is included, while in 
the right frame a matter- antimatter asymmetry A/A\> avy ~ 1/3 is assumed. Notice the 
much larger regions of acceptable parameter space in cases with asymmetry. We use the 
following parameters: Mi=300 GeV, M 2 =300 GeV, fi=Q00 GeV, tan (3 = 50 and m h =llh 
GeV. The region above the solid line is excluded by measurements of the invisible Z decay 
width jH]. 

in the expansion ov = a + bv 2 + 0(v A ), and thus is very efficient for annihilations at low 
velocities which take place in the late universe. At low velocities, the process vv — > uu, 
via t-channel neutralino exchange, is the most important annihilation channel. 

Given these characteristics, searches for gamma-rays, anti-protons, positrons or other 
particles produced in sneutrino annihilations in the Galactic halo are quite inefficient. 
Due to the rather large elastic scattering cross sections of sneutrinos, however, indirect 
detection techniques involving the capture of dark matter in bodies such as the Sun can 
be very effective. 

Neutrino telescopes detect neutrinos produced in dark matter annihilations in the 
core of the Sun by observing 'tracks' of muons generated in charged current interactions. 
If our sneutrino dark matter candidate is of muon or tau flavor, the flux of muons gener- 
ated in this way in a detector at Earth is on the order of 10 -12 cm" 2 s -1 jS], in conflict 
with data from the Super-Kamiokande and MACRO experiments, for example |17j . If the 
sneutrino candidate, however, is of electron flavor, these bounds can be easily averted. 

6 Conclusions 

In the standard freeze-out calculation for a weakly interacting dark matter relic, there 
is little reason to expect a density of dark matter which is similar to the density of 
baryons, short of anthropic arguments. In this letter, however, we have presented an 
explanation for the similarity between these two quantities. Our solution introduces 
an asymmetry between dark matter particles and anti-particles which is related to the 
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Figure 2: The thermal relic density as a function of mass for sneutrinos+antisneutrinos 
with no asymmetry (dot-dash), with a matter-antimatter asymmetry of A/A^ry — 1/6 
(solid) and the estimate of Eq.fjHJ) (dots). The relic density range favored by WMAP is 
bound by dashed lines (0.129 > n CDM /i 2 > 0.095). sin 0=0.3, Mi=300 GeV, M 2 =300 
GeV, /i=600 GeV, tan/? = 50 and m h =115 GeV have been used. See text for details. 
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baryon-antibaryon asymmetry. This leads to a natural dark matter relic density of the 
same order of magnitude as the baryon density. 

As an example, we considered a mixed sneutrino dark matter candidate which trans- 
fers its particle-antiparticle asymmetry to the baryons through the electroweak anomaly. 
We carry out the relic density calculation for such a candidate and find substantial re- 
gions of parameter space in which the observations of WMAP can be satisfied. With 
no asymmetry, only a narrow strip of parameter space can provide the observed relic 
density. 
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Abstract 

The inferred values of the cosmological baryon and dark matter densities are strikingly 
similar, but in most theories of the early universe there is no true explanation of this 
fact; in particular, the baryon asymmetry and thus density depends upon unknown, and 
a priori unknown and possibly small, CP-violating phases which are independent of all 
parameters determining the dark matter density. We consider models of dark matter 
possessing a particle-antiparticle asymmetry where this asymmetry determines both the 
baryon asymmetry and strongly effects the dark matter density, thus naturally linking J\ 
and fidm- We show that sneutrinos can play the role of such dark matter in a previously 
studied variant of the MSSM in which the light neutrino masses result from higher- 
dimensional supersymmetry-breaking terms. 



1 Introduction 



For some time it has been apparent that the inferred values of the cosmological baryon 
and dark matter densities are strikingly similar. For example, the latest WMAP- 
determined range for the dark matter density, 0.129 > ^cdm^ 2 > 0.095, is within a 
factor of a few of the combined WMAP and big-bang nucleosynthesis determined value 
of the baryon density, 0.025 > Q h h 2 > 0.012 [HI2]- I n the vast majority of models of 
the early universe, the cosmological baryon and dark matter densities are independently 
determined. The surviving baryon density is set by a baryon asymmetry generated dur- 
ing baryogenesis, and thus depends upon unknown, and a priori unknown and possibly 
small CP-violating phases, as well as unknown baryon-number violating dynamics. In 
contrast, the dark matter density is set by the 'freeze-out' of the interactions that keep 
the dark matter in equilibrium, and is independent of the dynamics of baryogenesis. In 
particular, although weakly-interacting massive particles at the TeV-scale naturally have 
a relic density of 0(1) times the critical density, this is not at all the case for baryons. 
The Boltzmann equations show the size of the baryon relic density in the absence of an 
asymmetry is O(10 -11 ). Thus the comparative closeness of the baryon and dark mat- 
ter densities poses a puzzle. One possible approach to this problem is to more closely 
integrate the dynamics of baryogenesis with that of the origin of dark matter. In partic- 
ular, it is natural to consider models where the dark matter and baryon sectors share a 
quantum number, either continuous or discrete, which provides a relation between their 
surviving number densities and thus energy densities. 1 

Specifically, in this letter we consider models of dark matter possessing a particle- 
antiparticle asymmetry where this asymmetry strongly effects the dark matter density, 
and through the electroweak (EW) anomaly, determines the baryon asymmetry, thus 
naturally linking and ^cdm- (A noteworthy early attempt along these lines which 
shares some features with our model is contained in Ref. jU].) The shared quantum 
number clearly requires that the dark matter particle not be it's own antiparticle, so the 
putative dark matter candidate cannot be the standard LSP neutralino of the broken 
Minimal Supersymmetric Standard Model (MSSM). In fact, we will show in sections 3 
and 4 that sneutrinos can play the role of such dark matter in a previously studied 
variant of the MSSM in which the light neutrino masses result from higher- dimensional 
supersymmetry-breaking terms ^UJEIEIEIIE]- This model preserves all the successes 
of the MSSM, such as stability of the weak scale and unification of gauge couplings, while 
being, at least in part, testable at the LHC. 

Before we proceed to the calculation of the relic dark matter and baryon densities 
and the details of our model, it is useful to consider some aspects of the idea of a shared 
quantum number determining the ratio of baryon to dark matter. First, consider the 
simplified case in which, to a very good approximation, the dark matter and baryon 
sectors cannot exchange their conserved quantum number after its first production. As 
an example, consider the situation in which the asymptotic baryon number B = 1 states 
of the Standard Model (SM) have global charge, q, while the lightest asymptotic states 

*An alternate approach is to invoke a form of the Anthropic Principle, see, for example, Ref. 
Other ideas are explored in Refs. [H □ El 
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(with mass m^m) in a hidden sector carry charge, Q. Then conservation of global charge 
implies that q(n h — n^) = —Q(rid m — n-^), where the n's are the number densities 
of the indicated (anti)particles. Further assume that interactions in both sectors are 
strong enough such that almost all antiparticles are eliminated by annihilation with 
their particles (assuming Q/q < for simplicity). This implies 

n b = cn dm with c=\Q/q\, (1) 

which in turn leads to 

7j— = = c • (2) 

Thus the energy densities are indeed related, but the ratio is only naturally 0(1) if the 
ratio of the baryon to dark matter masses is not too small. It is important to note that 
this disfavours models of the above type where the dark matter candidate is a 'hidden 
sector' particle and favours a particular class of models where the dark matter candidate 
has weak scale mass, its mass arising either from electroweak symmetry breaking, or from 
the dynamics which drives electroweak breaking, such as supersymmetry softly broken 
at the weak-scale. 2 

As a consequence, the naive statement of Eq.(^) also requires modification as the 
approximation of negligible interactions between the two sectors, and thus negligible ex- 
change of the shared conserved quantum number, does not hold in realistic models. In 
this situation, the ratio of the conserved quantum number stored in the two sectors, c, 
is not simply given by Q/q, but is determined by the 'chemical' equilibrium conditions 
between the two sectors just before the freeze-out of the relevant interactions. We will 
further discuss this point in section 2. Moreover, two other alterations to the naive re- 
lation Eq.(J2J) can be present. First, the shared quantum number may only be conserved 
modulo n, in which case 'self annihilations may occur, and second, it is not always the 
case that dark matter interactions are efficient enough such that n-^ <C fidm- Neverthe- 
less, as we argue in the next section, a careful consideration of such models shows that 
a version of the naive statement Eq.(J2J) can indeed hold. 



2 The Effect of a Matter- Antimatter Asymmetry on 
the Density of a Thermal Relic 

If there were an asymmetry between the density of dark matter particles and dark anti- 
matter particles in the early universe, it may have a considerable effect on the relic 
density calculation ^3]. In this section, we review the relic density calculation for a 
single species, two species with no asymmetry and two species with an asymmetry. 

If dark matter consists of a single species, the relic density calculation is straight- 
forward. The most common example of a single species dark matter candidate would be 

2 In extensions of the SM such as the MSSM where many gauge-non-neutral fermionic and scalar 
states also gain mass at this scale, the beta-function coefficient for SU(3)qct> sharply increases as 
energies are reduced below the weak scale, and the confinement scale and thus baryon mass is naturally 
not very much smaller than m wea k. 
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a neutralino LSP with negligible coannihilations with other SUSY particles. In such a 
case, the resulting relic density is simply 



, 1.04 x 10 9 x F 



VLh 2 ^ — ~ r, (3) 

Mpi ^gZ(a + 3b/x F y v ; 

where M-p\ is the Planck mass, xp is the WIMP's mass over the freeze-out temperature, 
g* is number of relativistic degrees of freedom available at freeze-out and a and b are 
related to the WIMP's self-annihilation cross section by ov = a + bv 2 + 0(v 4 ), v being 
the relative velocity of the annihilating WIMPs. For a weakly interacting relic, xp is 
generally about 20 and varies only mildly with changes in the cross section or mass. 

If we instead consider a dark matter candidate and its antiparticle, the calculation 
becomes somewhat different. Labeling these species 1 and 2, we now have each particle's 
self-annihilation cross section, an = a~22, and the particle-antiparticle annihilation cross 
section, 012 = 021 ■ In the limit that o\i is much smaller than an, we simply get two 
decoupled species which each provide the relic density of Eq.©. Otherwise, the effective 
cross section becomes an + 012 for each species. For the two species together, a density 
of 

„ t2 n 1.04 x 10 9 x F 1 

Qh 2 ^2 x - (4) 

Mpi v ^r(a 1 i + ai a + 3(6u + 6i 2 )/zf) 

is produced. This case is essentially the calculation of the density of a relic with the 
coannihilation of a degenerate particle [T51 . 

If an asymmetry exists between the number density of particles and antiparticles in 
the early universe, again this calculation becomes modified. The most well known exam- 
ple of this is the baryon-antibaryon asymmetry which leads to the current baryon density 
of the universe despite the very large baryon-antibaryon annihilation cross section. In 
such a case, where the particle-particle and antiparticle-antiparticle annihilation cross 
sections (011 and 022) are zero, or negligibly small, the minimum relic density can be 
related to the magnitude of the particle-antiparticle asymmetry: 

= ^Lv^ — — , (5) 



^bary 



where Qh1 &j:y is the current density of baryons in the universe, A and Ab ary are the 
particle-antiparticle asymmetries of our relic and of baryons, defined by A = (n — n)/n. 
^bary is known to be of the order of 10~ 9 . m and mb ary are the masses of our relic and 
of baryons (i.e. the proton mass). Of course, if the particle-particle annihilation cross 
section for the relic is not negligible, Eq. (0) will not hold. There will be a generic tendency 
for the density of a relic to move towards this value as a result of an asymmetry, however. 

To accurately assess the effect of an asymmetry on the thermal relic abundance, the 
effective cross sections must be determined by integrating over the thermal history. In 
particular, for each species separately, the densities predicted are 

Q ^ 2 ^ 1.04 x 10 9 x F 1 

Mpi Jgl (aii >e ff + ai2, e ff + 3(&u ieff + bi2 jeS )/x F ) 
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and 




„ , o 1.04 x 10 9 x F 1 

^ 2 = 77 F=7 ; 7^777 -r r~/ r > (7) 



-^Pl \/#* (a22,eff + a21,eff + 3(6 2 2,eff + b 2 l,^[)/xp 

where the effective cross sections are given by 

a ll,eff = All j a 12,eff = %fJ^^T 

022,eff = a 22, a 21,eff = XpJ^^T 

&ll,eff = &11> &12,eff — l 

W = & 22 , &2Mf = 24/~^^±4Jete. 

The factor, e~ x /(e~ x + A), is the depletion of the second species due to the asymmetry 
at a temperature, T = m/x. The effect of this depletion depends on the size of the 
asymmetry and the cross sections involved. If the estimate of Eq. © is considerably less 
than the standard two-species result of Eq.(J3J), then the asymmetry has little impact 
and the sum of Eqs.(JB} and (jZj) is nearly equal to Eq.(j3J. If, however, the estimate of 
Eq.fjSJ) is larger than the standard two-species result, the relic density moves from the 
value found with Eq.Q towards that of Eq.@. If crii/cr 12 is small, then the result will 
be closer to Eq.fjSJ). If (J\\l<J\i is large, then the density will be less modified. Note that 
in some cases, the effect of an asymmetry can be to lower the density of a relic by up to 
a factor of two. This occurs when o\\ and an are comparable thus efficiently deplete the 
second species while leaving the first largely uneffected. 

Although these effects are rather model dependent, the net effect can be large regions 
of parameter space which predict the observed density of dark matter for a WIMP with 
a mass near m ~ m p fiwMAph 2 /flbar y h 2 , for a matter-antimatter asymmetry similar to 
that for baryons. For concreteness, we consider a specific example of such a model in the 
following section. 



3 An Example: Mixed Sneutrino Dark Matter 

Within the context of the Minimal Supersymmetric Standard Model (MSSM), sneutrinos 
do not make a very appealing dark matter candidate. They annihilate quite efficiently, 
resulting in a relic density smaller than the observed dark matter density for sneutrinos 
lighter than about 500 GeV. Furthermore, their elastic scattering cross section is suffi- 
ciently large to be easily observed by direct dark matter experiments. Thus it is natural 
that in contrast to other supersymmetric dark matter candidates, sneutrino dark matter 
has received considerably less attention [T7j . 

Looking beyond the vanilla MSSM, however, it is quite easy to construct a viable 
model for sneutrino dark matter. For example, sneutrinos could mix through /L-terms 
with a singlet forming a mixed sneutrino state. This would reduce the annihilation cross 
section, potentially providing the appropriate quantity of dark matter, as well as reduce 
the elastic scattering cross section to experimentally acceptable levels. In fact, such a 
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scenario has previously been discussed as a way to reconcile the results of the direct dark 
matter detection experiments DAMA and CDMS [T2] • These models have light sneutrino 
mass eigenstates which are mixtures of left-handed 'active' sneutrinos, and right-handed 
'sterile' sneutrinos. 

Such mixing naturally arises in models in which the light neutrino masses are gener- 
ated by supersymmetry breaking effects fUl El GUI E] • The importance to us of these 
models is that the light sneutrino states share a non-anomalous (B — L)-symmetry with 
the baryons which is only weakly broken by the Majorana neutrino masses. In analogy to 
the Giudice-Masiero mechanism for the Higgs //-term there exists a global symmetry that 
forbids a direct mass, MrNN, and neutrino Yukawa coupling, XLNH U , involving the SM 
singlet right-handed neutrino superfield N, but allows non-renormalizable couplings to 
a supersymmetry-breaking hidden-sector field T. The right-handed neutrino mass and 
the neutrino Yukawa coupling, XLNH U , then only arise from l/M planck -suppressed terms 
involving the fundamental supersymmetry breaking scale mi ~ (vM p i) 1 / 2 , where v is the 
weak scale. The relevant terms in the Lagrangian are 

C = J d4e [ h M NN + h M^ N N + hB ^T NN + ■■■)+] d2 ° {^LNH U ) , (9) 

where all dimensionless couplings g, h, etc, are taken to be 0(1). The structure of the 
Lagrangian can be enforced by some simple symmetries [TTH ITT] IT3]. including R-parity, 
which leads to a stable lightest supersymmetric particle (LSP) dark matter candidate. 
After supersymmetry is broken in the hidden sector at the scale my, the F and A com- 
ponents of the field T acquire intermediate scale vevs, (Ft) ~ (At) 2 ~ m 2 . Substituting 
these vevs into the Lagrangian shows that after SUSY breaking, there arises: (1) a neu- 
trino Yukawa coupling of size ~ 1CT 7 — 1CT 8 , (2) right-handed neutrino and sneutrino 
masses at the weak scale, (3) a weak-scale trilinear scalar A-term, and (4) a right-handed 
sneutrino lepton-number violating i?-term with mag nitude B 2 ~ (few x 100 MeV) 2 . 
These together lead to two sources of neutrino masses, a tree level see-saw contribution, 
as well as a dominant 1-loop contribution, leading to a consistent and attractive neutrino 
phenomenology [TOJ E3] . 

More relevant for us, the sneutrino mass matrix squared is 

\ ( m2 T -h=Avsm(3 \ / p \ 
'* n) , * a ^ , ^ f 1 , (10) 



-j^Av sin p m\ 

where we have ignored the small lepton number violating I?-term for the moment. There- 
fore mixing between the left-handed 'active' sneutrino, u, and the right-handed 'ster- 
ile' sneutrino, n, occurs, producing the light mass eigenstate Pi = —vsinO + n* cos9, 
where tan 29 = — \/2Av sin/3/(m|. — m|), with m\ = + m 2 Rsoh and m\ = 
m| cos 2/3/2 + m 2 Lsoft . Since the coupling of the lighter sneutrino eigenstate to the Z 
is suppressed by sin#, the direct LEP experimental constraints are weakened, and the 
calculation of the sneutrino relic abundances is modified. A numerical analysis of the 
susy parameter space shows that a mixing sin 9 ~ 0.3 is quite natural. 3 (The B 2 mass- 
term splits the CP-even and CP-odd components of v\ by the amount 5 ~ cos 2 9B 2 jm\. 

3 Note that the size of the A terms is limited by considerations of vacuum stability |14|. 
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This can have important consequences for detection, potentially making consistent the 
DAMA and CDMS experiments depending upon parameters ^2]- For the purposes of 
our discussion this splitting is irrelevant.) 

We now turn to the calculation of the relative asymmetry in the sneutrino and 
baryon sectors. The method is a simple adaption of the standard 'chemical' equilibriation 
techniques applied to, for example, the calculation of the ratio B/(B — L) in the SM [19 
or MSSM j20] i n the presence of anomaly-induced baryon number violating processes in 
the early universe. For each species of particle we introduce a chemical potential, [Aa, 
and for each interaction in equilibrium, say A + B <-> C + D there is a corresponding 
relation among chemical potentials \la + A*b — A*c + A*d- Note that when in thermal 
(not necessarily chemical) equilibrium the CP-conjugate species A has ii-^ = —\ia- The 
asymmetry in particle-antiparticle number density for the relevant case of \x <C T and a 
species of mass m is given by n — n = gT 2 fi/6 for T ^> m, and is Boltzmann suppressed 
by exp(— m/T) for T < m. 

The electroweak B + L anomaly interactions provide an efficient mechanism to trans- 
mute between B- and L-number densities, at least above the EW phase transition tem- 
perature T c where they are rapid. Our task is to calculate the ratio of sneutrino to baryon 
asymmetries at T c where the exchange between the two sectors switches off. (We implic- 
itly assume that the EWPT is 1st order so that anomaly-induced interactions drop out 
of equilibrium below T c . For realistic values of T c the results of our calculation are not 
particularly sensitive to this assumption.) Above the EWPT temperature the chemical 
potentials for all SM gauge bosons vanish, so, for example, members of an EW doublet 
have the same \i. Moreover, efficient flavour mixing among quarks induced by EW inter- 
actions such as uu + 0i <-> URji involving the neutral up- like higgs, or du + 4>2 ^Rj-> 
involving the neutral down-like higgs, ensures that only one chemical potential, say that 
for ul, is required to describe the quark sector. Similarly, efficient flavour mixing among 
the light neutrinos holds for mixing angles, 8ij > 10~ 7 , so that only one chemical poten- 
tial, (j,„ L , is needed to describe the neutrino and sectors. In addition non-perturbative 
EW processes can be thought of creating, say, a neutron made of lhd quarks, and a lhd 
neutrino out of the vacuum for each generation. When this process is in equilibrium it 
imposes an additional relation which can be written 

3fi UL + fi UL = 0. (11) 

When the MSSM susy spectrum is included the analysis becomes more involved. To 
simplify the analysis let us first assume at a temperature T (with T > T c ) that the 
MSSM susy spectrum, including k rhd sneutrinos can be considered light (m < T). As 
shown in Ref. [20] the usual soft-susy-breaking terms and supersymmetric partners of the 
SM interactions lead to enough relations among the chemical potentials that only 3 are 
needed for an independent set: liq = (/^ 10 — /^ 20 )/2 in the neutral Higgs sector, lm Ul , and 
fi UL . In our case we must add an additional k light (m < T c ) rhd neutrinos and sneutrinos 
with their associated potentials \i VRi and fiv m , but we also have A-term interactions 
between the lhd and rhd sneutrinos, and the very small lepton-number violating rhd 
sneutrino B-term. These together enforce the conditions 
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where the first equality arises from A-term mediated interactions (as the A terms are 
0{m w ) it is simple to show that these are in equilibrium), and the second equality arises 
from SU(2) W x U(l)y gaugino interactions. 4 In terms of //o,£Vl, an d Hv L , the condition 
Eq.()ll|l is unchanged in the case of the MSSM + k rhd neutrino supermultiplets. 

The final step in the calculation of the relative number densities proceeds by ex- 
pressing the relevant densities in terms of the independent chemical potentials. For the 
baryon number density B the result is simply 

B = 6T 2 /v- (13) 

For the sneutrino asymmetry the result depends upon the mass spectrum of the sneutrino 
mass eigenstates relative to T c . In the range of parameters of interest to us where the 
light sneutrino eigenstate u% is the LSP we have mi < T c . On the other hand the heavier 
eigenstate 1)2 can have m 2 > T c or m 2 < T c . If both eigenstates satisfy m < T c then the 
sneutrino density is (allowing for a possible k such states) 

kT 2 2kT 2 

S = —(^i>L + »Ur) = —^-VVL ( 14 ) 

while if only mi < T c then S is half as large. Finally using the non-perturbative relation 
Eq.JIU we find tnat 

t (1^) 

^bary B 3 6 

depending on m 2 . We view k = 1 as most likely, so we specialize to that case in the next 
section. Finally we mention that a full susy-mass-threshold dependent calculation of the 
relative asymmetries introduces only small corrections to this result (see Ref . |2Uj for a 
similar analysis in the MSSM case). 

In the end the vital point is that independent of whatever dynamics at scales E > T c 
produces an asymmetry in either the baryon sector, or in the neutrino or sneutrino 
sector, the (B + L)-anomaly-induced interactions together with EW gaugino and A-term 
interactions automatically distributes the asymmetry between the baryons and the light 
sneutrino states, with a predictable A/A^j ratio. The asymmetry could originate from, 
e.g. a GUT-based baryogenesis mechanism, or maybe more interestingly in the context of 
the sneutrino dark matter model there is the possibility that it could arise via a resonant 
leptogenesis mechanism at the TeV-scale as discussed in Ref.|14j. The end result is that 
we expect 1/3 > A/A^y > 1/6 independent of the source of the asymmetry. 



4 Asymmetric sneutrino and baryon abundances 

We will now apply the formalism of section El to the sneutrino model. Given the value of 
1/3 > A/Abary ^ 1/6, we expect the role of the matter-antimatter asymmetry to be most 

4 Note that the rhd neutrinos only interact via the small neutrino Yukawa, A ~ 10~ 7 — 10~ 8 , so 
these states are not in equilibrium at T ~ T c . However, fi l/R does not enter our final expression for the 
sneutrino to baryon number density ratio, so its value is irrelevant. 
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important for sneutrinos with masses on the order of m ~ few x m p fiwMAph 2 /flbaryh 2 — 
30 GeV. For this reason, we will focus on the mass range of about 10 to 80 GeV. In this 
range, the most important annihilation channel is s-channel Z-exchange to fermions. 
This is a p-wave amplitude, and has no a-term in the expansion in v 2 (see section |2J). 
This is also only a sneutrino-antisneutrino process (i.e. crn^, not cxi2,2i)- Other channels 
include t-channel neutralino and chargino exchanges to neutrinos and charged leptons. 
The channel to charged leptons includes both a and b terms in the v 2 expansion. Addi- 
tionally, sneutrino-sneutrino pairs (or antisneutrino-antisneutrino pairs) can annihilate 
via t-channel neutralino exchange to neutrino pairs of the same flavor. This is the only 
contribution to 012,21 in this model. 

The annihilation cross sections for the neutralino and chargino exchange processes 
depend on the masses and mixings of these exchanged particles, although the Z-exchange 
channels do not. We have adopted the parameters M!=300 GeV, M 2 =300 GeV, /i=600 
GeV and tan (3 = 50 throughout our calculation. Additionally, there is the mixing angle 
between the sneutrino and the singlet which can be varied. This quantity has the effect 
of reducing all electroweak couplings of the sneutrino eigenstate by a factor of sin 2 9, i.e. 
reducing the annihilation cross sections by a factor of sin 4 9. 

Finally, one can consider the s-channel exchange of Higgs bosons to bb (or ZZ and 
W + W~ above the appropriate mass thresholds). This process is not suppressed by sin 4 9. 
The efficiency of this process depends on the magnitude of the trilinear scalar coupling 
and is only important near the Higgs resonance (m„ ~ m^/2) or above the gauge boson 
mass thresholds. 

Our results are shown in figure H The shaded regions of the parameter space predict 
a relic density within the range measured by WMAP (0.129 > ^cdm^ 2 > 0.095). In 
the left frame, no matter-antimatter asymmetry was included. This frame can be easily 
understood as requiring greater suppression of electroweak couplings (smaller sin 9) as the 
sneutrino mass approaches the Z mass. The dip at 56-59 GeV is due to s-channel higgs 
exchange to bb. In the right frame, a matter-antimatter asymmetry of A/A^j ~ 1/6 was 
included. To further illustrate this effect, we show the result of this calculation across 
one value of sin 9 in figure El Below about 30 GeV, the matter-antimatter asymmetry 
has little effect on the calculation and the solid and dot-dashed lines fall nearly on top 
of each other. In the range of roughly 30-70 GeV, however, the asymmetry pulls the relic 
density above the standard symmetric result into the range favored by WMAP. Above 
this range, sneutrino-antisneutrino annihilation decreases, leading to larger relic densities 
for the case with no asymmetry. The relic density for the asymmetric case, however, is 
largely determined by the sneutrino-sneutrino annihilation cross section in this region, 
so does not increase as rapidly, therefore resulting in a relic density much closer to the 
preferred value, even for > 70 GeV. 

5 Indirect Detection 

Sneutrinos in the mass range we have studied here annihilate dominantly through s- 
channel Z exchange in the early universe. This annihilation channel consists of a 6-term 
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Figure 1: The regions of parameter space which provide the quantity of mixed sneutrino 
cold dark matter measured by WMAP, 0.129 > fi C DM^ 2 > 0.095. In the left frame, the 
standard calculation with no matter-antimatter asymmetry is used. In the center frame 
the effect of a matter-antimatter asymmetry with A/A^ ary ~ 1/6 is included, while in 
the right frame a matter- antimatter asymmetry A/A\> avy ~ 1/3 is assumed. Notice the 
much larger regions of acceptable parameter space in cases with asymmetry. We use the 
following parameters: Mi=300 GeV, M 2 =300 GeV, fi=Q00 GeV, tan (3 = 50 and m h =llh 
GeV. The region above the solid line is excluded by measurements of the invisible Z decay 
width [PS]. 

in the expansion av = a + bv 2 + O (t> 4 ) , and thus is very inefficient for annihilations at low 
velocities which take place in the late universe. At low velocities, the process vv — > uu, 
via t-channel neutralino exchange, is the most important annihilation channel. 

Given these characteristics, searches for gamma-rays, anti-protons, positrons or other 
particles produced in sneutrino annihilations in the Galactic halo are quite inefficient. 
Due to the rather large elastic scattering cross sections of sneutrinos, however, indirect 
detection techniques involving the capture of dark matter in bodies such as the Sun can 
be very effective. 

Neutrino telescopes detect neutrinos produced in dark matter annihilations in the 
core of the Sun by observing 'tracks' of muons generated in charged current interactions. 
If our sneutirno 5 dark matter candidate is of muon or tau flavor, the flux of muons gen- 
erated in this way in a detector at Earth is of the order of 10~ 12 cm -2 s _1 [T2j, in conflict 
with data from the Super-Kamiokande and MACRO experiments, for example [21]. If 
the candidate, however, is of electron flavor, these bounds can be easily averted. 

6 Conclusions 

In the standard freeze-out calculation for a weakly interacting dark matter relic, there 
is little reason to expect a density of dark matter which is similar to the density of 
baryons, short of anthropic arguments. In this letter, however, we have presented an 

5 In the model presented, the dark matter candidate is actually an antisneutrino and consequently 
antineutrinos will be produced by dark matter annihilations in the sun. 
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Figure 2: The thermal relic density as a function of mass for sneutrinos+antisneutrinos 
with no asymmetry (dot-dash), with a matter-antimatter asymmetry of A/A^ry — 1/6 
(solid) and the estimate of Eq.fjHJ) (dots). The relic density range favored by WMAP is 
bound by dashed lines (0.129 > n CDM /i 2 > 0.095). sin 0=0.3, Mi=300 GeV, M 2 =300 
GeV, /i=600 GeV, tan/? = 50 and m h =115 GeV have been used. See text for details. 
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explanation for the similarity between these two quantities. Our solution introduces 
an asymmetry between dark matter particles and anti-particles which is related to the 
baryon-antibaryon asymmetry. This leads to a natural dark matter relic density of the 
same order of magnitude as the baryon density. 

As an example, we considered a mixed sneutrino dark matter candidate which trans- 
fers its particle-antiparticle asymmetry to the baryons through the electroweak anomaly. 
We carry out the relic density calculation for such a candidate and find substantial re- 
gions of parameter space in which the observations of WMAP can be satisfied. With 
no asymmetry, only a narrow strip of parameter space can provide the observed relic 
density. 



Acknowledgments 



We wish to thank Yuval Grossman for a useful correspondence. SW is supported by 
PPARC Studentship Award PPA/S/S/2002/03530. 



References 



[1] D. N. Spergel et ai, Astrophys. J. Suppl. 148 (2003) 175 |arXiv:astro- ph/0302209| 



[2] S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B 592, 1 (2004). 
[3] F. Wilczek, arXiv:hep-ph/0408167| 

[4] A. Kusenko and M . E. Shaposhnikov, Phys. Lett. B 418, 46 (1998) 
|arXiv:hep-ph/9709492| . 



10 



[5] K. Enqvist and J . McDonald, Nucl. Phys. B 538, 321 (1999) 
|arXiv:hep-ph/9803380| . 

[6] D. H. Oaknin and A. Zhitnitsky, |arXiv:hep-ph/0309086| 

[7] D. H. Oaknin and A. R. Zhitnitsky, [arXiv:hep-ph/0406146| 

G. R. Farrar and G. Zaharijas, arXiv:hep-ph/0406281 

D. B. Kaplan, Phys. Rev. Lett. 68, 741 (1992). 



[9 
[10 

[11 
[12 

[13 
[14 

[15 

[16 
[17 



[18 
[19 
[20 



N. Arkani-Hamed, et al, Phys. R ev. D 64, 1150 11 (2001) |a rXiv:hep-ph /0006312|; 
N. Arkani-Hamed, et al, |arXiv:hep-ph/000700lj 

F. Borzumati and Y. Nomura, Phy s. Rev. D 64, 53005 (2001) 
|arXiv:hep-ph/0007018|; F. Borzumati, et al, |arXiv:hep-ph/0012118l 



D. R. Smith and N. Weiner, Phys. Rev. D 64, 043502 (2001) |arXiv:hep-ph/0101138 
D. Tucker-Smith and N. Weiner, arXiv:hep-ph/0402065, 

J. March-Russell and S. M. West, ar Xiv:hep-ph/0403067| 

T. Hambye, J. M arch-Russell an d S. M. West, JHEP 0407 (2004) 070 
|arXiv:hep-ph/0403183| ; S. M. West, |arXiv:hep-ph/04083"l8l 



K. Griest and D. Seckel, Nucl. Phys. B 283, 681 (1987) [Erratum-ibid. B 296, 1034 
(1988)]. 

K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

T. Falk, K. A. Oliv e and M. Srednicki, Phys. Lett. B 339, 248 (1994) 
|arXiv:hep-ph/9409270| ; L. E. Ibanez, Phys. Lett. B 137, 160 (1984); J. S. Hagelin, 

G. L. Kane and S. Raby, Nucl. Phys. B 241, 638 (1984); L. J. Hall, T . Moroi and 

H. Murayama, Phys. Lett. B 424, 305 (1998) larXiv:hep-ph/9712515] ; T. Han and 
R. Hempfling, Phys. Lett. B 415, 161 (1997) |arXiv:hep-p~/9708264T 



T. Hebbeker, Phys. Lett. B 470, 259 (1999) arXiv:hep-ph/9910326|. 

J. A. Harvey and M. S. Turner, Phys. Rev. D 42, 3344 (1990). 

T. Inui, T. Ichihara, Y . Mimura and N. Sakai, Phys. Lett. B 325, 392 (1994) 
|arXiv:hep-ph/9310268| . 



[21] G. Giacomelli [MACRO Collaboration], Braz. J. P hys. 33, 211 (2 003) 
|arXiv:hep-ex/0210006| ; I. De Mitri [MACRO collaboration] , |arXiv: hep-ex/ 02 12055) 



11 



